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2-D  SIMULATIONS  OF  THE  NRL  LASER  EXPERIMENT 


I.  Introduction 

We  present  the  results  of  two  sets  of  simulations  of  the  DNA/NRL  laser  experiment 
using  a  two-dimensional  gas-dynamic  code.  The  set  of  simulations  discussed  were 
prompted  by  the  structure  seen  in  the  expanding  shock  front  in  recent  experiments 
( e.g .  Stamper  et  al .,  1984  ).  These  experiments  were  run  in  the  high  ambient 
pressure,  hydrodynamic  regime.  In  general  these  experiments  show  the  formation 
of  expanding  bubbles  or  aneurysms  on  the  shock  front.  Data  from  the  experiment 
show  that  the  shock  has  the  character  of  a  self-similar  blast  wave;  such  a  blast  wave 
is  generally  stable.  The  question,  then,  that  these  simulations  attempt  to  answer 
is  whether  it  is  possible  to  perturb  the  initial  conditions  sufficiently  in  the  ablated 
target  material  or  in  the  ambient  surroundings  to  give  rise  to  the  aneurysms  seen  in 
the  experiment.  In  the  first  set  of  simulations  we  perturb  or  deform  the  distribution 
of  the  ablated  target  debris  in  order  to  discover  whether  the  resulting  blast  wave 
becomes  unstable  to  an  aneurysm-like  deformation.  In  the  second  set  of  simulations, 
we  study  how  an  inhomogeneity  in  the  ambient  initial  conditions  -  in  this  case  a 
pre-heated  channel  -  might  affect  the  propagation  of  the  blast  wave. 

A  second  question  can  be  studied  with  these  simulations  that  is  of  interest  for 
the  seeding  of  late-time  structure  in  a  HANE.  It  is  predicted  theoretically  (  e.g., 
Richtmeyer,  1960  )  and  observed  experimentally  (  Andronov  et  al.,  1976  )  that  the 
contact  surface  between  two  media  can  be  made  Rayleigh- Taylor  unstable  by  the 
passage  of  a  shock  across  the  boundary.  Such  a  process  can  occur  in  the  debris-air 
interface  of  a  HANE  and  might  be  responsible  for  phenomena  seen  at  later  times. 
Thus,  in  these  simulations  we  were  also  interested  in  the  question  of  whether  the 
applied  perturbation  in  the  debris  shell  damped  out  as  the  shocked  region  grew,  or 
whether  the  perturbation  grew  as  a  function  of  time. 

Manuscript  approved  October  17, 1984. 


n.  Perturbed  Debris  Shells 


In  this  section,  the  ambient  material  into  which  the  debris  moves  is  assumed  to 
be  uniform,  and  it  is  the  debris  shell  which  is  given  perturbations  of  various  forms 
to  see  whether  any  of  these  perturbations  grow  in  time.  Two  different  types  of 
perturbations  are  discussed  below.  The  first  type  is  what  is  usually  meant  by  a 
perturbation,  a  small  change  in  a  uniform  state.  Here  both  the  shell  position  and 
mass  distribution  were  given  such  small  variations.  The  second  type  of  perturba¬ 
tion  follows  a  suggestion  from  C.  Longmire  that  a  slower  jet  of  material  coming 
behind  the  main  debris  shell  may  be  responsible  for  the  aneurysms.  The  results 
of  simulating  such  jets  are  given  below.  We  use  the  term  perturbation  for  these 
simulations  loosely,  since  the  change  in  the  initial  conditions  here  is  quite  large  in 
magnitude. 

All  the  calculations  discussed  in  this  section  were  made  using  a  cylindrical 
geometry.  A  spherical  geometry  is  more  in  keeping  with  the  general  symmetry  of 
the  laser  experiment.  To  test  the  general  stability  of  debris  shells  we  prefer  not  to 
choose  a  special  coordinate  direction  as  we  would  if  we  used  the  polar  angle  in  a 
spherical  coordinate  system.  The  other  angular  coordinate  in  the  spherical  case  is 
the  azimuthal  one.  This  leads  to  reasonable  results  only  if  the  polar  angle  subtended 
is  quite  small  and  there  is  no  structure  in  that  direction.  Neither  of  these  conditions 
appears  to  be  well  satisfied  for  the  laser  experiment.  Thus,  cylindrical  symmetry 
was  chosen  as  a  compromise  that  could  be  used  over  large  angular  extents  and  still 
contain  the  effects  of  the  radial  divergence  terms  in  the  equations.  The  following 
set  of  equations  were  solved: 
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where 


E  =  2p(v°*  +  +  -  l 

Pe  =  pv$ 

Pr  =  PVr 

-y  =  ratio  of  specific  heats 


The  equations  were  solved  in  conservation  form  using  an  extension  of  the  par¬ 
tial  donor  cell  method,  PDM  (  Hain  1978,  Lyon,  1984).  Eighth  order  centered 
spatial  differencing  and  a  second  order  time  integration  scheme  (modified  Euler) 
were  employed.  The  numerical  grid  was  40  x  40  cells. 

The  grid  was  continuously  adapted  in  the  r  direction  during  the  calculation, 
while  the  6  grid  was  fixed.  The  endpoints  of  the  mesh  in  the  radial  direction  were 
moved  to  keep  all  structure  on  the  grid.  The  interior  mesh  points  were  spaced  to 
provide  resolution  in  the  regions  of  steep  gradients.  This  was  done  by  minimizing 
the  integral, 
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at  each  time  step.  The  radial  coordinate,  r,  is  assumed  to  be  a  function  of  the 
mesh  coordinate,  *.  The  second  term  in  the  integrand  is  a  smoothness  measure  - 
the  smaller  the  integral  the  smoother  the  mapping  of  r  to  ».  The  first  term  weights 
i  to  regions  where  w  is  large.  In  these  calculations  the  weight  function, 


was  used.  This  attempts  to  concentrate  points  in  the  regions  of  large  density  gra¬ 
dients.  The  quantity,  A,  is  an  adjustable  parameter  used  to  provide  control  over 
the  balance  between  the  two  terms.  For  the  series  of  runs  discussed  here  A  was 
adjusted  so  that  a  maximum  grid  compression  about  three  times  that  of  uniform 
spacing  occurred  in  the  debris  shell. 

The  initial  conditions  were  designed  to  mimic  shot  number  13261  in  the  NRL 
laser  experiment.  This  particular  shot  shows  a  single  clearly  defined  aneurysm  ( 
Stamper  et  al.,  1984).  The  ambient  gas  pressure  was  5  Torr  with  an  ambient 
mean  molecular  weight  about  equal  to  that  of  the  aluminum  target.  The  debris 
velocity  and  mass  ablated  for  the  shot  energy  of  38  J  were  estimated  from  the 
scaling  relations  given  by  Ripin  (1983).  The  debris  velocity  used  was  575  km/sec. 
The  shell  width  was  taken  to  be  4  ns  times  575  km/sec,  or  .23  cm.  The  shell  was 
initially  centered  at  0.2  cm  from  the  target,  and  it  was  assumed  that  there  was 
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no  ambient  material  within  the  shell’s  maximum  radius.  The  mass  ablated  was 
estimated  as  .23  fx g,  which  gave  an  estimated  shell  density  of  7  x  1O-0  gm/cm3. 
This  is  an  estimate  of  the  shell  density  for  the  spherical  case,  but  was  used  for  these 
cylindrical  simulations. 

A.  Perturbed  Shells 

Two  types  of  perturbations  were  applied  to  the  initial  debris  shell.  In  the  first,  the 
position  of  the  shell  was  given  a  small  sinusoidal  ripple  in  the  6  direction  amounting 
to  about  5%  of  the  thickness  of  the  shell.  In  the  second,  the  shell  position  was 
unperturbed  but  the  mass  density  within  the  shell  was  given  a  sinusoidal  perturba¬ 
tion.  In  all  cases  no  aneurysms  were  observed  to  form  and  any  disturbance  of  the 
position  of  the  outgoing  shock  quickly  damped  out. 

The  case  of  a  10%  sinusoidal  mass  perturbation  is  typical.  Figure  1  shows  the 
density  contours  in  the  simualtion  at  two  different  times.  In  the  first,  Figure  l.a, 
2.58  ns  after  the  beginning  of  the  calculation,  we  can  see  a  double  shock  structure 
where  the  debris  shell  has  interacted  with  the  ambient  gas.  One  shock  is  propagating 
outward  into  the  ambient  gas,  the  other  inward  through  the  debris  shell.  There  is 
also  a  contact  discontinuity  separating  the  ambient  from  debris  material.  It  is  not 
apparent  in  the  density  plot  because,  by  coincidence,  the  initial  shell  density  and 
ambient  density  were  the  same.  Figure  l.b  shows  the  density  about  3  ns  later.  The 
outer  shock  shows  a  slight  bulge  in  the  region  of  enhanced  density  -  the  center  of 
the  plot.  This  bulge  quickly  damps  out,  however,  as  seen  in  Figure  l.c.  This  figure 
shows  the  density  some  6.5  ns  later.  The  blast  wave  front  now  appears  to  be  very 
uniform;  the  initial  bulge  has  disappeared.  In  all  the  runs  performed  with  these 
types  of  simulations  the  blast  wave  front  was  very  stable  to  anything  resembling  an 
aneurysm. 

While  these  runs  gave  stable  blast  waves  there  was  a  considerable  amount  of 
internal  structure  in  the  region  of  the  contact  surface  between  the  debris  material 
and  the  ambient  gas.  The  developing  irregularities  can  be  seen  in  Figure  1,  but 
it  is  easier  to  see  them  in  Figure  2  which  gives  a  contour  plot  of  the  temperature. 
Figure  2.a  shows  the  temperature  at  about  6  ns  into  the  calculation.  The  two  shock 
fronts  are  apparent,  and  also  apparent  is  a  rippling  of  the  contact  surface  region 
in  between  the  two.  Figure  2,b  gives  the  temperature  6  ns  later.  A  very  large  and 
growing  distortion  can  be  seen  now  in  this  same  region.  The  rate  of  growth  appears 
to  be  consistent  with  Richtmeyer’s  (  1960  )  theory  which  predicts  linear  rather  than 
exponential  growth.  Thus,  a  cylindrical  or  spherical  geometry  probably  does  not 
severely  affect  this  mechanism.  Irregularities  in  the  debris  shell  can  cause  large 
irregularities  in  the  internal  structure  of  the  blast  wave  bubble. 


B.  Slow  Jets 

These  simulations  were  intended  to  test  the  idea  that  aneurysms  could  be  formed  by 
a  jet  of  material  leaving  the  target  at  a  slower  speed  than  the  rest  of  the  ejected  ma¬ 
terial.  When  the  shock  wave  produced  by  the  initial  material  has  slowed  sufficiently, 
the  slower  material  can  overtake  it  and  go  on  ahead  producing  an  aneurysm-like 
protuberance  of  the  blast  wave  surface.  In  these  rims,  the  middle  10  cells  of  the 
calculation  in  the  6  direction  were  given  an  enhanced  density  and  a  decreased  ve¬ 
locity.  The  densities  and  velocities  were  chosen  so  that  the  “jet”  region  always 
had  a  higher  momentum  and  energy  than  the  rest  of  the  shell.  With  this  partic¬ 
ular  set  of  initial  conditions,  the  bulk  of  the  shell  should  move  more  quickly  than 
the  “jet”  until  enough  mass  is  swept  up  for  the  blast  wave  solution  to  be  set  up. 
Thereafter,  the  higher  energy  in  the  “jet”  should  allow  it  to  overtake  the  rest  of  the 
shell  and  proceed  ahead.  If  there  were  no  interaction  between  the  two  regions  of 
the  shell,  at  late  times  both  would  follow  a  blast  wave  solution  with  differing  pro¬ 
portionality  constants.  This  would  imply  that  if  this  mechanism  were  responsible 
for  aneurysms,  the  aneurysms  would  be  self-similar  structures  on  the  surface  of  the 
blast  wave.  However,  there  clearly  is  strong  coupling  between  the  two  regions,  and 
it  is  the  purpose  of  these  calculations  to  determine  whether  that  interaction  helps 
or  hinders  the  running  ahead  of  the  “jet” . 

Figure  3  shows  the  density  evolution  of  a  “small”  jet.  In  this  run,  the  density 
enhancement  was  6p  jp  =  0.4  and  the  velocity  decrease  was  Sv  jv  =  -0.1.  The 
increase  in  jet  energy  was  thus  «  20%.  In  Figure  3.a,  2.6  ns  after  the  start  of 
the  calculation,  the  enhanced  density  jet  can  be  seen  lagging  slightly  behind  the 
the  rest  of  the  shell.  In  Figure  3.b,  3  ns  later,  the  shock  in  front  of  the  increased 
density  region  has  caught  up  with  the  rest  of  the  shell.  The  increased  density  is 
still  apparent  in  the  contours  at  this  time.  Figure  3.c  shows  the  density  at  a  very 
much  later  time,  39  ns.  There  is  a  slight  irregularity  in  the  shock  front,  but  the 
jet  region  now  lags  behind  the  rest  of  the  shell.  Also  the  density  of  the  shell  has 
become  much  more  uniform.  It  is  difficult  to  pick  out  an  increased  density  region  in 
the  6  driection,  except  very  far  inward  from  the  shock.  The  transition  from  Figure 
3.b  to  3.c  is  a  smooth  one.  In  intermediate  times  the  shock  front  in  the  jet  region 
does  run  slightly  ahead  of  the  rest  of  the  shell,  but  the  major  evolution  seems  to  be 
a  decrease  in  the  density  non-uniformity  behind  the  shock  in  the  9  dizection. 

A  higher  density  jet  was  also  simulated  as  shown  in  Figure  4.  Here  6p  /p  = 
1.0  and  6v  /v  =  -0.1.  Thus,  the  energy  within  the  shell  is  almost  twice  that  of  the 
rest  of  the  shell.  On  the  assumption  that  both  regions  follow  their  own  blast  wave 
solutions,  this  jet  should  run  ahead  of  the  rest  of  tha  shell  by  about  15%  of  the 
shock  radius.  Figure  4. a  shows  the  density  very  early  in  the  calculation.  The  jet 
is  already  forging  slightly  ahead  of  the  rest  of  the  debris.  In  figures  4,b  and  4.c  at 
times  6.1  and  13.4  ns,  respectively,  the  process  continues.  Note,  however,  that  the 


density  irregularity  behind  the  shock  is  strongly  decreasing.  Figure  4.c  is  about  the 
peak  of  the  advance  of  the  jet  over  the  rest  of  the  material.  The  final  figure  4.d 
shows  the  density  at  41.4  ns.  Here,  the  jet  area  lags  behind  the  rest  of  the  shell. 
Note  also  the  density  clumps  behind  the  shock.  The  strange  behavior  of  the  flanks 
running  ahead  of  the  jet  is  almost  certainly  due  to  the  periodic  boundary  conditions 
used  in  the  d  direction.  As  the  jet  moves  on  ahead  of  the  rest  of  the  shell  the  shock 
wave  in  front  of  it  slows  relative  to  the  rest  of  the  shell  and  also  broadens.  As  this 
braodening  reaches  the  boundaries,  a  cusp-like  structure  is  formed  as  the  shock  front 
and  its  image  meet  at  an  oblique  angle  at  the  boundary.  This  dumps  shock  energy 
in  the  region  of  the  boundary,  and  the  shock  in  the  boundary  region  accelerates 
and  overtakes  the  jet  region.  While  such  an  analysis  is  not  strictly  valid  for  an 
expanding  cylindrical  shock,  an  explanation  in  terms  of  shock  geometric  optics  can 
be  found  in  Whitham  (  1974  ).  Once  the  boundaries  run  ahead,  they  cause  the 
same  process  to  occur  in  the  center  of  the  grid  ,  so  that  an  oscillatory  behavior 
is  set  up.  This  oscillation  eventually  damps.  While  interesting,  the  phenomenon 
does  depend  on  the  artificial  boundary  conditions  used,  and  so  is  probably  of  little 
practical  interest. 

The  density  clumps  observed  behind  the  shock  may  be  of  considerable  interest 
for  a  HANE,  however.  We  note  that  tney  begin  to  form  even  at  t  =  6.1  ns,  far 
too  early  for  the  boundary  conditions  to  affect  the  simulation.  At  later  times  the 
disturbance  appears  to  cascade  to  a  relatively  short  angular  wavelength.  In  the  final 
plot,  Figure  4.d,  the  amplitude  of  these  clumps  is  about  2  to  1  over  the  surroundings. 
The  mechanism  producing  these  clumps  is  not  known,  although  the  initial  clump 
of  to  the  side  of  the  jet  in  Figure  4.b  is  apparently  related  to  the  vorticity  set  up  by 
the  differing  velocities  in  the  jet  and  the  rest  of  the  shell.  It  is  quite  possible  that 
the  large  amplitude  to  which  these  clumps  grow  is  due  to  the  periodic  boundary 
conditions.  However,  these  small  scale,  high-amplitude  irregularites  could  be  a 
significant  source  of  structure  in  at  least  some  HANE’s. 


HI.  The  Pre-Heated  Channel 

From  the  previous  section  it  is  apparent  that  an  energy  conserving  explosion  in  the 
fluid  regime  is  not  spontaneously  unstable  to  aneurysm-like  modes.  In  fact,  very 
special  initial  conditions  with  large  amplitude  irregularities  must  be  formulated  to 
produce  anything  resembling  aneurysms.  In  this  section  we  treat  an  alternative 
hypothesis;  We  assume  that  something  creates  a  channel  with  a  narrow  angular 
extent  in  front  of  the  blast  wave.  The  conditions  for  the  propagation  of  the  shock 
within  the  channel  are  then  different  from  those  outside.  The  shock  could  then  run 
faster  within  the  channel  than  outside,  leading  to  an  aneurysm-like  appearance.  A 
weak  piece  of  evidence  in  favor  of  this  idea  is  that  in  the  current  series  of  laser  shots, 
the  aneurysm  is  usually  single  and  roughly  aligned  with  the  laser  or  target  normal 
axis.  In  fact,  in  one  shot,  where  the  target  was  tilted  relative  to  the  laser  beam, 
two  aneurysms  were  seen  -  one  along  the  laser  axis  and  the  other  along  the  target 
normal.  We  can  only  speculate  about  the  possible  cause  of  channel  formation,  but 
two  possibilities  are  a  “hot  spot”  in  the  laser  beam  and  a  beam  of  high  energy  ions 
produced  by  some  microphysical  process  in  the  laser-target  interaction. 

There  are  a  number  of  ways  in  which  a  pre-heated  channel  can  affect  the 
propagation  of  a  blast  wave.  One  possibility  is  that  the  shock  wave  runs  faster 
through  hotter  material.  Unfortunately,  though  the  effect  is  in  the  proper  direction 
to  produce  aneurysms,  the  magnitude  is  rather  small.  The  speed  of  strong  shocks  is 
to  lowest  order  independent  of  the  temperature  of  the  ambient  gas.  The  first  order 
correction  to  that  speed,  6vjv,  is  «  7  —  2//yM2,  where  M  is  the  Mach  number  of 
the  shock  in  the  ambient  gas.  The  smaller  the  Mach  number,  the  faster  the  shock 
runs  for  7  <  2.  However,  even  if  the  gas  is  pre-heated  to  a  temperature  of  about 
one  eV,  the  corection  still  would  appear  to  be  only  about  1%.  A  more  sizeable  effect 
is  required.  A  second  possibility  is  that  if  the  gas  is  not  in  local  thermodynamic 
equilibrium  (  LTE  ) ,  the  effective  *y  of  the  gas  could  change,  thus  affecting  the  shock 
speed.  There  is  a  relatively  strong  dependence  of  shock  speed  on  ^  -  Rah  oc  ~i°‘95 
for  7  =  1.2.  While  this  is  an  interesting  possibility,  simulation  of  this  case  requires 
a  very  complex  numerical  model  including  chemical  and  transfer  effects.  Thus,  we 
do  not  pursue  this  possibility  here. 

A  third  possibility  is  the  one  upon  which  we  will  concentrate  in  this  section. 
A  pre-heated  channel  can  give  rise  to  a  lowered  density  ahead  of  the  blast  wave. 
The  process  is  described  as  follows:  Suppose  a  cylindrical  channel  is  heated  at  some 
instant  near  the  time  of  the  laser  pulse.  There  is  a  significant  amount  of  time 
between  that  instant  and  the  instant  when  the  blast  wave  overruns  the  channel. 
During  this  time  the  heated  gas  is  undergoing  expansion  away  from  the  channel  axis. 
If  the  material  within  the  channel  is  much  hotter  than  its  surroundings,  the  density 
in  the  channel  can  become  much  lower  than  ambient.  The  rough  time  scale  for  such 
expansion  is  tc  =  /?rc/cc,  where  (3  is  a  constant  on  the  order  of  unity,  and  rc  and  cc 


are  the  initial  radius  and  sound  speed  of  the  channel,  respectively.  For  the  ambient 
gas  in  ther  laser  experiments,  assuming  rc  =  0.5mm  and  Tc  =  5  eV  gives  tc  =  (3  (100 
ns).  The  time  scale,  tc  is  thus  roughly  in  the  right  range  to  be  significant  in  the  laser 
experiment.  The  quantity,  0,  is  greater  than  one  for  a  channel  in  which  there  is  a 
step-function  radial  dependence  of  heating.  The  rarefaction  wave  which  evacuates 
the  center  does  not  arrive  there  until  tc.  A  Gaussian-like  radial  dependence  leads 
to  faster  results.  Numerical  experiments  have  shown  that  a  significant  lowering  of 
density  (  «  20%  )  at  the  center  of  the  channel  occurs  at  about  0.3  tc.  In  this  case 
tc  is  defined  by  ac/cc(r  =  0),  where  ac  is  the  dispersion  of  the  Gaussian  and  cc 
is  evaluated  at  thje  center  of  the  channel.  A  halving  of  the  central  density  occurs 
at  about  2tc.  The  density  dependence  of  the  blast  speed  upon  density  is  not  very 
strong.  Still,  assuming  that  the  propagation  of  the  shock  down  the  channel  is  given 
by  its  own  spherical  blast  wave  dependence,  t.e. 
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we  can  see  that  a  20%  density  reduction  leads  to  about  a  4%  speed  increase.  That  is 
not  a  large  enough  difference  to  produce  aneurysms  by  itself,  but  it  is  large  enough 
to  make  non-linear  effects  important. 

Before  going  on  to  the  simulations,  it  is  worthwhile  to  discuss  one  further  back- 
of-the-envelope  argument.  The  formation  of  aneurysms  occurs  over  a  wide  range  of 
laser  energies  and  ambient  densitiies.  As  a  rule  of  thumb,  the  lower  the  density  or 
the  higher  the  laser  energy,  the  later  does  the  aneurysm  formation  occur.  We  can 
see  whether  the  pre-heated  channel  mechanism  is  consistent  with  such  scaling  by 
setting  tc  equal  to  the  time  it  takes  the  blast  wave  to  reach  a  given  distance  from 
the  origin.  That  is, 
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where  the  radius  of  the  channel,  rc,  is  &cRah  and  k  is  a  proportionality  constant.  As 
noted  above  the  shock  radius  is  proportional  to  f0,4.  The  radius  of  the  channel  scales 
at  a  power  closer  to  one.  For  very  long  times  we  might  expect  the  channel  expansion 
to  go  over  into  a  cylindrical  blast  wave,  where  rc  would  scale  as  f0,5.  In  any  event, 
the  implication  is  that  for  all  radii  greater  than  the  one  found  by  solving  (5),  the 
channel  should  be  well  formed.  Since  we  are  only  looking  for  scaling  arguments,  we 
solve  (5)  for  a  constant  channel  energy  dropping  most  of  the  constants  to  obtain 
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?his  shows  a  fairly  weak  dependence  on  laser  energy  and  ambient  density,  but 
he  dependencies  are  in  the  right  direction  -  higher  laser  energy  and  lower  density 
ead  to  slower  aneurysm  formation.  It  might  be  reasonable  to  assume  that  the 
>re- heating  mechanism  has  an  inverse  square  dependence  from  the  target  spot. 
Assuming  Ech  =  E0Rq/R 2,  we  get 
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which  has  a  linear  dependence  on  E\,i  and  pQl.  Once  again,  however,  the  depen¬ 
dencies  axe  in  the  right  direction  to  be  consistent  with  the  laser  results.  These 
arguments  support  the  idea  that  the  pre-heated  channel  is  an  explanation  of  the 
aneurysm  phenomenon.  The  next  sections  deal  with  a  numerical  simulation  of  the 
effects  of  a  pre-heated  channel. 


A.  Simulations 


A  slightly  different  form  of  the  code  used  in  §  II  was  used  to  simulate  the  pre-heated 
channel.  While  still  two-dimensional,  the  variables  were  radius,  r,  from  t1  e  target 
and  the  polar  angle,  0.  This  is  the  most  reasonable  coordinate  sytem  to  use  for  a 
2-D  simulation  as  long  as  there  is  a  preferred  axis  along  which  the  perturbation  is 
expected  to  grow.  Such  is  the  case  for  the  prc-hcated  channel.  Secondly,  while  a 
cylindrically  symmetric  blast  wave  does  not  behave  too  differently  from  a  spherical 
one,  a  slab  symmetric  channel  (which  would  be  the  case  for  a  cylindrical  blast  ) 
acts  quite  differently  from  a  cylindrically  symmetric  one.  Since  there  seems  to  be 
a  fairly  sensitive  trade-off  between  time  scales  associated  with  the  blast  wave  and 
the  channel,  it  is  important  in  this  case  to  use  an  accurate  geometry. 

The  equations  solved  for  these  simulations  were 
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Fig.  6  —  Density  contours  at  selected  times  for  a  simulation  with  a  constant  initial  channel  temperature. 
The  display  format  is  similar  to  Figure  1.  The  contours  are  labeled  in  units  of  1  x  10“ 7  gm/cm3.  The 
times  displayed  are:  6,a  *  64.1  ns,  and  6.b  =  89.1  ns. 


Fig.  5  (Cont’d)  —  Pressure  contours  at  selected  times  for  a  simulation  with  a  constant  initial  channel 
temperature.  The  display  format  is  similar  to  Figure  1.  The  contours  are  labeled  in  units  of  1  x  106 
ergs/cm3.  The  times  displayed  are:  5.a  —  48.9  ns,  5.b  *■*  64.1  ns,  and  5.c  **  89.1  ns. 
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(c) 


DENSITY  T=4 1 .39 


(d) 

Fig.  4  (Cont’d)  —  Density  contours  at  selected  times  for  a  "large"  jet  type  of  perturbation.  Display 
format  is  similar  to  Figure  1.  Contours  are  labeled  in  units  of  10~7  gm/cm3.  The  times  displayed  are: 
4,a  »  2.6  ns,  4.b  —  6.09  ns,  4.c  —  13.39  ns,  and  4.d  *»  41.4  ns. 
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DENSITY  T=  2.62 
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(b) 


Fig.  4  —  Density  contours  at  selected  times  for  a  "large"  jet  type  of  perturbation.  Display  format  is 
similar  to  Figure  1.  Contours  are  labeled  in  units  of  10" 7  gm/cm3.  The  times  displayed  are:  4.a  —  2.6 
ns,  4.b  —  6.09  ns,  4.c  *  13.39  ns,  and  4,d  *  41.4  ns. 
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Fig.  3  (Cont’d)  —  Density  contours  at  selected  times  for  a  "small"  jet  type  of  perturbation.  Display 
format  is  similar  to  Figure  1.  Contours  are  labeled  in  units  of  10~ 7  gm/cm3.  The  times  displayed  are: 
3. a  —  2.6  ns,  3.b  —  5.89  ns,  and  3.c  —  39.3  ns. 
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Fig.  2  —  Similar  to  Figure  1  except  the  temperature  is  displayed.  The  times 
displayed  are:  2. a  ™  5  83  ns  and  2.b  —  1 1.89  ns. 
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Fig.  1  (Cont’d)  —  Density  contours  at  selected  times  for  a  debris  shell  with  an  initial  sinusoidal  density 
perturbation.  The  data  is  contoured  on  an  angular  wedge  shaped  area.  The  bottom  axis  gives  the  radial 
distance  along  that  edge  of  the  wedge.  Contours  are  labeled  in  units  of  1(T7  gm/cm3.  The  times 
displayed  are:  l.a  —  2.58  ns,  l.b  —  5.86  ns,  l.c  —  12.27  ns. 
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Fig.  1  —  Density  contours  at  selected  times  for  a  debris  shell  with  an  initial  sinusoidal  density 
perturbation.  The  data  is  contoured  on  an  angular  wedge  shaped  area.  The  bottom  axis  gives  the  radial 
distance  along  that  edge  of  the  wedge.  Contours  are  labeled  in  units  of  10-7  gm/cm3.  The  times 
displayed  are:  l.a  —  2.58  ns,  l.b  —  5.86  ns,  l.c  —  12.27  ns. 


Giuliani  (1984)  has  suggested  an  interesting  mechanism  for  the 
aneurysms  which  depends  on  the  thermal  conductivity  being  high  in  the  region  of 
the  aneurysm  and  low  elsewhere.  The  most  plausible  reason  for  this  is  a  magnetic 
field  which  is  present  everywhere  in  the  shell  except  the  aneurysm.  Keskinen  (1983) 
has  pointed  out  how  self-generated  magnetic  fields  can  be  produced  in  the  initial 
laser  expansion.  The  same  mechanism  can  also  operate  within  the  region  that 
bounds  the  main  shell  from  the  lower  density  region  produced  by  the  channel.  The 
rate  at  which  self-generated  magnetic  fields  are  produced  is  ( e.g .  Tidman,  1975) 
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This  requires  that  the  pressure  gradient  and  the  density  gradient  be  in  somewhat 
different  directions.  Reference  to  Figures  5  and  6  show  that  in  the  channel  region 
it  is  quite  possible  for  the  gradients  to  be  misaligned.  The  misalignment  can  be  as 
high  as  15-20°.  Typical  gradient  scale  lengths  for  the  pressure  and  density  are  of 
the  order  of  .05  cm  or  less  at  the  channel  boundary.  With  p  «  2  x  10~5  gm/cm 
and  P  as  108  erg/cm3,  we  get  dB/dt  ~  10 12  G/sec  for  a  10°  misalignment.  Taken 
at  face  value,  this  result  would  imply  that  kiloGauss  fields  could  be  generated  at 
the  channel  boundary.  This  is  quite  likely  a  gross  over-estimate.  However,  it  does 
indicate  that  self-generated  magnetic  fields  may  be  quite  important  for  the  evolution 
of  the  pre-heated  channels.  In  particular,  it  gives  a  natural  way  in  which  Giuliani’s 
mechanism  can  be  seeded  and  maintained,  since  the  self-generated  magnetic  field 
at  the  channel  boundary  will  suppress  thermal  conduction  from  the  channel  region 
to  the  rest  of  the  shell. 
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ambient  medium.  The  temperature  at  the  center  of  the  channel  was  given  by 


2 

Tch  =  200  eV  )  exp(— (r  —  0.3  cm)/0.3  cm) 

\0.3  cm/ 

This  gives  a  very  high  energy  of  200  eV  at  the  edge  of  the  debris  region  initially. 
The  energy  drops  very  rapidly,  however,  to  less  than  2  eV  at  1  cm  from  the  target. 
The  intent  here  was  to  see  if  a  channel  could  be  opened  at  very  early  times  in 
the  calculation,  and  once  opened  whether  the  bulge  created  could  be  maintained 
or  grow  as  the  shock  ran  into  less  evacuated  regions.  Figure  7  shows  the  pressure 
profiles  at  49,  67,  and  95  ns  into  the  calculation.  Figure  7.a  (  49  ns  )  shows  a  bulge 
that  is  consistent  with  the  size  of  aneurysms  seen  at  that  stage  of  the  blast  wave 
expansion.  On  the  other  hand,  the  succeeding  two  figures  show  the  bulge  damping 
while  simulataneously  broadening.  It  is  clear  that  once  the  perturbation  in  ambient 
density  into  which  the  blast  wave  runs  dies  out,  so  does  the  bulge  in  the  shock  front. 


IV.  Discussion 

Two-dimensional  gas-dynamic  simulations  of  the  NRL  laser  experiment  have  been 
performed  to  study  the  formation  of  aneurysms  in  the  blast  wave,  and  to  study  the 
formation  of  structure  internal  to  the  blast  wave  itself.  Two  sets  of  simulations  have 
been  performed:  perturbed  debris  shell  and  pre-heated  channel. 

The  conclusions  from  the  perturbed  debris  shell  simulations  are  summarized 
as  follows.  The  source  of  the  aneurysms  observed  in  the  NRL  laser  experiment  is 
probably  not  due  to  perturbations  of  the  debris  shell.  No  “simple”  perturbation 
of  the  shell  leads  to  the  appearance  of  a  simple  bump  growing  on  the  surface  of 
the  blast  wave.  The  heuristic  arguments  about  the  stability  of  blast  waves  seem 
to  be  vindicated  by  this  set  of  simulations.  On  the  other  hand,  a  large  amount 
of  structure  can  be  produced  relatively  easily  internal  to  the  shock  front.  With 
suitable  perturbations  it  appears  that  density  enhancements  of  a  factor  of  two  may 
be  possible.  Further  work  is  needed,  however,  to  determine  how  general  such  large 
density  jumps  may  be. 

On  the  other  hand,  the  simulations  using  a  pre-heated  channel  show  that 
such  a  channel  can  effectively  produce  a  distortion  in  the  spherical  shape  of  an 
expanding  blast  wave.  However,  it  appejjs  that  the  magnitude  of  the  distortion  is 
too  small  to  explain  the  aneurysms.  Another  problem  with  the  pre-heated  channel 
explanation  is  that  the  aneurysms  seem  to  have  the  highest  density  of  any  region  of 
the  shock  front.  The  simulations  show  the  bulges  having  somewhat  lower  density 
than  the  surrounding  shock.  It  is  unlikely  that  a  pre-heated  channel  can  be  the  full 
explanation  of  the  aneurysms. 
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B.  Results 


In  the  first  simulation  the  channel  was  assumed  to  have  a  constant  angular  width 
and  the  same  pressure  everywhere  along  the  central  axis  initially.  The  pressure 
corresponded  to  a  thermal  energy  of  4  eV  per  particle.  This  is  perhaps  the  largest 
energy  that  could  be  deposited  in  a  channel  without  producing  much  higher  values 
than  are  seen  by  the  spectroscopic  diagnostics  which  sample  a  long  line  of  sight 
through  which  the  channel  might  run.  This  particular  type  of  initial  condition 
favors  the  development  of  the  aneurysm  at  late  times  in  the  calculation.  This  is 
born  out  by  Figure  5,  which  shows  the  time  development  of  the  pressure  in  this 
calculation.  We  plot  the  pressure  because  it  gives  a  better  indication  of  the  shock 
strength  than  does  the  density.  In  Figures  5.a  -  5.c  at  times  49,  64,  and  89  ns  into 
the  calculation,  respectively,  we  can  see  a  slight  protuberance  of  the  shock  front 
developing  and  increasing  in  both  in  distance  ahead  of  the  rest  of  the  shock,  and, 
most  particularly,  in  its  angular  extent.  Figure  6,  which  plots  the  density  at  64 
and  89  ns,  shows  the  same  effect.  Also  apparent  is  the  decrease  in  density  in  the 
channel.  In  Figure  6.b  the  minimum  density  in  the  channel  just  ahead  of  the  shock 
is  almost  a  factor  of  2  below  the  ambient.  Even  at  the  earlier  time,  the  density 
decrease  is  about  40%.  It  is  clear  that  the  position  of  the  shock  front  in  the  channel 
does  not  show  the  increase  in  speed  that  would  be  expected  if  were  following  its  own 
blast  wave  shock  radius-time  relation.  Interestingly  enough,  however,  if  we  solve 
the  blast  wave  relation  for  the  ratio  of  P/p  as  a  function  of  time,  we  find  that  it 
goes  as 


«  Vsh  * 


(9) 


The  jump  in  the  density  ratio  from  channel  to  ambient  is  2.  The  pressure  ratio 
jump  just  behind  the  shock  is  about  5/7.  This  is  quite  close  to  the  value  that  would 
be  predicted  by  (9)  by  changing  the  density  by  a  factor  of  two.  In  this  particular 
respect,  then,  the  shock  propagation  within  the  channel  does  follow  the  blast  wave 
relation.  The  reason  that  the  shock  within  the  channel  does  not  run  way  ahead  of 
the  the  rest  of  the  channel  is  that  the  bulge  has  a  significantly  larger  surface  area 
than  the  area  on  a  sphere  with  the  unperturbed  blast  wave  radius  which  subtends 
the  same  solid  angle  from  the  center  of  the  explosion.  The  energy  of  the  shock  is 
thus  spread  over  a  larger  area,  and  the  shock  speed  slows  down  over  what  it  would 
be  in  a  strictly  spherical  expansion.  The  fact  that  the  pressure  jump  is  about  right 
for  a  blast  wave  solution  suggests  that  the  volume  into  which  the  set  amount  of 
energy  flows  is  about  the  same  as  it  would  be  for  a  purely  spherical  expansion. 

The  second  simulation  in  this  sequence  usqd  a  different  initialization  of  the 
channel  to  attempt  to  mimic  a  diverging  beam  which  was  also  absorbed  by  the 
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E  =  ±p(ve2  +  tv2)  + 

P6>  =  pve 
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*7  =  ratio  of  specific  heats 

and  the  polar  angle,  9,  is  taken  to  be  the  angle  from  the  axis  of  the  channel.  The 
same  rezoning  technique  used  in  §  II  was  used  again  here  to  increase  resolution  in 
the  region  of  the  shock  front. 

We  will  discuss  two  simulations  below,  both  of  which  assume  a  laser  energy  of 
10  J  and  an  ambient  pressure  of  5  Torr.  The  initial  conditions  used  have  for  R  <0.3 
cm  from  the  target 

P  —Pa  +  Pd  =  9  x  10“ 6  gm/cm3 

vr  = — — — Vd  =  128  km/sec 
Pa  +  Pd 

P  .fcLLjjgyi  ( i  -  -W-)  =  5  x  10°  erg/cm3 
2  V  Pa+PdJ  ' 

where  the  a  and  d,  subscripts  refer  to  ambient  and  debris,  respectively.  The  debris 
velocity  was  assumed  to  be  the  575  km/sec  used  in  the  previous  section  and  the 
debris  mass  was  scaled  down  by  the  ratio  of  10/38,  giving  a  debris  density  «2x 
10“b  gm/cm  .  It  is  assumed  that  ambient  and  debris  material  are  mixed  uniformly 
and  have  fully  coupled.  For  R  >  0.3  cm  from  the  target,  the  ambient  gas  is  assumed 
to  have  a  temperature  of  0.5  eV  except  in  the  channel.  The  energy  deposited  in  the 
channel  is  assumed  to  have  a  Gaussian  6  dependence  -  Ech  oc  exp  (—62/29^)  -  and 
in  both  simulations  9C  =  0.07  radians.  In  the  first  simulation  a  constant  channel 
energy  of  4  eV  was  assumed,  while  in  the  second  the  form  used  atttempted  to  mimic 
a  spreading  beam  with  absorption  by  the  ambient  gas. 
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Fig.  7  (Cont’d)  —  Pressure  contours  at  selected  times  for  a  simulation  with  an  initialization  to  mimic  a 
diverging  beam  with  ambient  absorption.  The  display  format  is  similar  to  Figure  1.  The  contours  are 
labeled  in  units  of  1  x  106  ergs/cm3.  The  times  displayed  are  7.a  —  49.4  ns,  7.b  —  66.5  ns,  and  7.c  — 
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ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

01 CY  ATTN  DICK  STEINHOF 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

01 CY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 

0 1 CY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA,  VA  22311 
01 CY  ATTN  J.M.  AEIN 
01 CY  ATTN  ERNEST  BAUER 
01  CY  ATTN  HANS  WOLFARD 
01  CY  ATTN  JOEL  BENGSTON 

INTL  TEL  &  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

01 CY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

11011  TORREYANA  ROAD 

P.O.  BOX  85154 

SAM  DIEGO,  CA  92138 

0 1 CY  ATTN  J.L.  SPERLING 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

01  CY  ATTN  DOCUMENT  LIBRARIAN 
01CY  ATTN  THOMAS  POTEMRA 
01  CY  ATTN  JOHN  DASSOULAS 


KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
0 1 C Y  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED 
STUDIES 

816  STATE  STREET  (P.O  DRAWER  30) 

SANTA  BARBARA,  CA  93102 
01 CY  ATTN  DASIAC 
01 CY  ATTN  WARREN  S.  KNAPP 
0 1 C  Y  ATTN  WILLIAM  MCNAMARA 
01 CY  ATTN  B.  GAMBILL 

LINKA8IT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 

0 1 CY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  8  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
01 CY  ATTN  DEPT  60-12 
0 1 C  Y  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  8  SPACE  CO.,  INC. 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

01 CY  ATTN  MARTIN  WALT  DEPT  52-12 
01 CY  ATTN  W.L.  IMHOF  DEPT  52-12 
01 CY  ATTN  RICHARD  G.  JOHNSON, 

DEPT  52-12 

01 CY  ATTN  J.B.  CLADIS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 

01 CY  ATTN  R.  HEFFNER 

M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 

01 CY  ATTN  DAVID  M.  TOWLE 
01 CY  ATTN  L.  LOUGHLIN 
01 CY  ATTN  D.  CLARK 


MCDONNEL  DOUGLAS  CORPORATION 
5301  80LSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
01  CY  ATTN  N.  HARRIS 
0 1 C Y  ATTN  J.  MOULE 
0 1 C Y  ATTN  GEORGE  MROZ 
0 1 C  Y  ATTN  W.  OLSON 
0 1  CY  ATTN  R.W.  HALPRIN 
0 1 C  Y  ATTN  TECHNICAL 

LIBRARY  SERVICES 

MISSION  RESEARCH  CORPORATION 

735  STATE  STREET 

SANTA  BARBARA,  CA  93101 


PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

01 CY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01801 

01 CY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 


01  CY 

ATTN 

P.  FISCHER 

BELLEVUE 

,  WA 

98009 

01  CY 

ATTN 

W.F.  CREVIER 

01  CY 

ATTN 

E.J.  FREMOUW 

01  CY 

ATTN 

STEVEN  L.  GUTSCHE 

01  CY 

ATTN 

R.  BOGUSCH 

PHYSICAL 

DYNAMICS,  INC. 

01  CY 

ATTN 

R.  HENDRICK 

P.O.  BOX 

10367 

01  CY 

ATTN 

RALPH  KIL8 

OAKLAND, 

CA 

94610 

01  CY 

ATTN 

DAVE  SOWLE 

ATTN  A 

.  THOMSON 

0 1  C  Y 

ATTN 

F.  FAJEN 

01  CY 

ATTN 

M.  SCHEIBE 

R  8  D  ASSOCIATES 

01  CY 

ATTN 

CONRAD  L.  LONGMIRE 

P.O.  BOX 

9695 

01  CY 

ATTN 

B.  WHITE 

MARINA  DEL  REY,  CA  90291 

01  CY 

ATTN 

R.  STAGAT 

01  CY 

ATTN 

FORREST  GILMORE 

01  C  Y 

ATTN 

WILLIAM  B.  WRIGHT, 

MISSION 

RESEARCH  CORP.  ' 

0 1  C  Y 

ATTN 

WILLIAM  J.  KARZAS 

1720  PANDOLPH 

1  ROAD,  S.E. 

01  CY 

ATTN 

H.  ORY 

ALBUQUERQUE, 

NEW  MEXICO  87106 

0 1  C  Y 

ATTN 

C.  MACDONALD 

01  CY 

R.  STELLINGWERF 

01CY 

ATTN 

R.  TURCO 

01  CY 

M.  ALME 

01  CY 

ATTN 

L.  DeRAND 

01  CY 

L.  WRIGHT 

01  C  Y 

ATTN 

W.  TSAI 

MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

0 1 C  Y  ATTN  JOHN  MORGANSTERN 
0 1 C  Y  ATTN  G.  HARDING 
01 CY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1320  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
01 CY  ATTN  W.  HALL 
01 C  Y  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
1 2  3  AO  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

0 1 C  Y  ATTN  E.C.  FIELD,  JR. 


RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
01 CY  ATTN  CULLEN  CRAIN 
0 1 CY  ATTN  ED  8EDR0ZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

01 CY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42nd  STREET 
NEW  YORK,  NY  10036 

0 1 C Y  ATTN  VINCE  TRAPANI 
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SCIENCE  APPLICATIONS,  INC. 
1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

01 CY  ATTN  LEWIS  M.  LINSON 
0 1 C  Y  ATTN  DANIEL  A.  HAMLIN 
0 1 C  Y  ATTN  E.  FRIEMAN 
0 1 C  Y  ATTN  E.A.  STRAKER 
01CY  ATTN  CURTIS  A.  SMITH 

SCIENCE  APPLICATIONS,  INC 
1710  GOODRIDGE  DR. 

MCLEAN,  VA  22102 
0 1 C  Y  J.  COCKAYNE 
01 C  Y  E.  HYMAN 


SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 


MENLO 

PARK, 

CA 

94025 

01  CY 

ATTN 

J  . 

CASPER 

0 1  C  Y 

ATTN 

DONALD  NEILSON 

01  CY 

ATTN 

ALAN  BURNS 

01  CY 

ATTN 

G. 

SMITH 

01  CY 

ATTN 

R. 

TSUNODA 

01  CY 

ATTN 

DAVID  A.  JOHNSON 

01  CY 

ATTN 

WALTER  G.  CHESNUT 

01  CY 

ATTN 

CHARLES  L.  RINO 

01  CY 

ATTN 

WALTER  JAYE 

01  CY 

ATTN 

J  . 

VICKREY 

01  CY 

ATTN 

RAY 

L.  LEADABRAND 

01  CY 

ATTN 

G. 

CARPENTER 

01  CY 

ATTN 

G. 

PRICE 

01  CY 

ATTN 

R. 

LIVINGSTON 

01  CY 

ATTN 

V. 

GONZALES 

01  CY 

ATTN 

D. 

MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

0 1 C  Y  ATTN  W.P.  BCQUIST 

TOYON  RESEARCH  CO. 

P.O.  Box  6890 
SANTA  BARBARA,  CA  93111 
01  CY  ATTN  JOHN  ISE,  JR. 

01 C Y  ATTN  JOEL  GAR8ARIN0 

TRW  DEFENSE  8  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90273 
01CY  ATTN  R.  K.  PLEBUCH 
01 CY  ATTN  S.  ALTSCHULER 
01CY  ATTN  D.  DEE 
01 CY  ATTN  D/  STOCKWELL 
SNTF/1575 


VISIDYNE 

SOUTH  BEDFORD  STREET 
8U  RUNG  TON,  MASS  01803 
01 CY  ATTN  W.  REIDY 
01 C  Y  ATTN  J.  CARPENTER 
0 1 C  Y  ATTN  C.  HUMPHREY 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 

01 CY  ATTN:  N.  ZABUSKY 


DIRECTOR  OF  RESEARCH 
U.S.  NAVAL  ACADEMY 
ANNAPOLIS,  MD 
02CY 


34 


END 

FILMED 

7-85 

DTIC 


